Andreev current induced by ferromagnetic resonance 
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We study charge transport through a metallic dot coupled to a superconducting and a ferromag- 
netic lead with a precessing magnetization due to ferromagnetic resonance. Using the quasiclassical 
theory, we find that the magnetization precession induces a dc current in the subgap regime even 
in the absence of a bias voltage. This effect is due to the rectification of the ac spin currents at the 
interface with the ferromagnet; it exists in the absence of spin current in the superconductor. When 
the dot is strongly coupled to the superconductor, we find a strong enhancement in a wide range of 
parameters as compared to the induced current in the normal state. 
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Spin-transfer torque allows one to manipulate the mag- 
netization of a ferromagnetic (F) layer by means of a 
spin-polarized current [l], 0. Random-access memories 
using this effect in order to induce magnetization rever- 
sal of the active elements are on their way to commer- 
cialization. The reverse effect, namely the generation of 
a spin current in a normal metal (N) by means of a dy- 
namically precessing ferromagnetic metal, has also been 
predicted [3|. In the absence of direct spin probes, this 
effect may be measured by using a second ferromagnet as 
an analyzer that converts the spin current into a charge 
current. However, it was pointed out theoretically Q and 
measured experimentally that a single F/N junction 
is enough to both generate and detect the spin current 
through the generation of a dc voltage at ferromagnetic 
resonance (FMR) in an open-circuit geometry. At the 
origin of this phenomenon is the spin accumulation on 
the normal side of the junction - due to the precession- 
induced spin current - which is typically different from 
that on the ferromagnetic side. If transmissions for the 
majority and minority electron species through the junc- 
tion are different, the difference in spin accumulation gen- 
erates a net charge current which must be compensated 
by a difference in electrochemical potentials such that 
no charge accumulation occurs. Spin relaxation inhibits 
spin accumulation and, thus, suppresses the effect. 

The aim of our work is to explore how this effect is 
modified in a ferromagnet / superconductor junction. The 
combination of ferromagnetic and superconducting (S) 
materials has been shown to lead to a variety of inter- 
esting spin phenomena @, [!]■ However, the study of the 
interplay between magnetization dynamics and supercon- 
ductivity is a relatively new topic. Experimentally, a nar- 
rowing of the FMR width at the superconducting transi- 
tion was observed in an F/S bilayer @. Theoretically, it 
was proposed that a dynamically precessing ferromagnet 
may generate a long-range proximity effect 110]. This 
effect would manifest itself in the enhancement of the 
critical current in a phase-biased ferromagnetic Joseph- 
son junction under FMR conditions. Related signatures 
in the tunneling density of states of the F layer have also 
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FIG. 1: Setup of the junction. A metallic dot is coupled 
to a ferromagnetic lead with precessing magnetization m(t) 
on the left and to a normal or superconducting lead on the 
right. The left barrier is characterized by the conductances 
Gi and G m , defined above Eq. Q, whereas the right barrier 
is characterized by the conductance G r - 



been investigated However, these works disregard 



interface effects and, therefore, do not take into account 
the possible generation of an FMR-induced dc voltage. 
Finally, let us note that the first experiments on voltage 
generation by FMR of Refs. @, i] were performed with 
Al as the normal metal, which becomes superconducting 
at low temperatures. 

As the FMR-generated charge current in an F/N junc- 
tion is typically associated with a spin current, one may 
wonder what happens in an F/S junction in the sub- 
gap re gim e, where transport is mediated by Andreev pro- 
cesses |l2j. We show that the generation of charge cur- 
rent in the absence of a spin current in a conventional 
singlet superconductor is possible. In fact, the absence 
of spin currents in the superconductor may even lead to 
a strong enhancement of the induced charge current as 
compared to the normal state. 

As in the normal state, the two main ingredients nec- 
essary to generate the effect are spin-dependent trans- 
missions through the junction and a spin accumulation 
region 19!) . The simplest setup meeting these require- 
ments is a metallic dot coupled through tunnel barriers 
to a ferromagnet and to a superconductor (see Fig. [T]) . 
Coulomb blockade effects are neglected, assuming that 
the conductances of the barriers largely exceed the con- 
ductance quantum. 

The magnetization precession in the ferromagnetic 
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lead is described by a time-dependent exchange field, 
J(t) = Jm(t) with 



The Green function in the F lead, gi , is determined by 



m(t) = (sin 9 cos fit, sin 9 sin fit, cos 9), 



(1) 



acting on the spin of the conduction electrons. Here the 
precession frequency, fi, and the tilt angle, 9, are both 
tunable with external dc and rf fields under standard 
FMR conditions [l3| . We will consider them as externally 
fixed parameters. 

The precession of the magnetization drives the system 
out of equilibrium and, thus, may generate a current. To 
describe the system, we use the quasiclassical Keldysh 
theory Q. In particular, the current through the junc- 
tion can be expressed in terms of the quasiclassical Green 
function g of the dot. Here g is a matrix in Keldysh, 
Nambu, and spin space. In Keldysh space, it has a trian- 
gular structure with retarded (g R ) 7 advanced (g A ), and 
Keldysh (g K ) components. Furthermore, it satisfies the 
normalization condition g 2 = 1. 

The equations determining the Green functions are 
most conveniently written in the rotational frame for 
the magnetization precession, where the problem is sta- 
tionary Assuming that the conductance of the dot 
largely exceeds the conductances of the junctions with 
the leads, the equation determining g may be cast in the 
form 
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+ Il+ I r = 0. (2) 



Here <7j and tj are Pauli matrices in Nambu and spin 
space, respectively (i = x,y,z). Furthermore, Gq — 
e 2 /ir is the conductance quantum (in units where H = 1), 
and S is the mean level spacing in the dot. The spin- 
dependent energy shift ±fi/2 is a spin-resolved chemical 
potential induced by the transformation from the labo- 
ratory to the rotational frame. The boundary conditions 
with the ferromagnetic (I = left) and superconducting 
(r = right) leads are represented by the matrix currents 
lu r and depend on the Green functions gu r describing 
the non-equilibrium state in the leads due to the magne- 
tization precession. 

Tunneling through an F/N interface is generally spin- 
dependent. The relevant processes can be characterized 
by the total conductance of the junction, Gi, and the 
difference between the conductances for the majority and 
minority electrons, G m [20]. The matrix current at the 
tunnel interface between the dot and the ferromagnet 



then takes the form [15 1 



It = y + -ril™ ■ <"~z,gi},!j], 



(3) 



where m = m(0). Within the quasiclassical approxima- 
tion, we assume \G m \ <C G;. Thus, G m can be treated 
perturbatively. 



fi 
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(4) 



where the self-energy E = ~iTg^(E + (fi/2)<7 z ) accounts 
for inelastic scattering in the relaxation time approxi- 
mation. Here, 1/T is the inelastic scattering time and 
<7at is the equilibrium Green function in a normal metal. 
Namely, g* (A) (E) = ±r z and g%(E) = 2tJ(E), where 
f(E) — tanh(£'/2T) is related to the Fermi distribution 
at temperature T. For a large exchange field, J 3> fi,T, 
the solution of Eq. ((4]) takes the form gf^ — ±r z 
and gf — 2r z (f + + /_ cos 9m ■ a), where f±(E) = 
[/(£ + fi/2)±/(£-fi/2)]/2. 

The matrix current at the dot-superconductor tunnel 
interface is given as 



G 1 



[gr,g], 



(5) 



where G r is the conductance of the junction. The Green 
function in the S lead reads g r — gs(E+ (fi/2)cr z ), where 
gs is the equilibrium Green function in a superconductor. 
Namely, gf A) {E) = (-iEr z + At x )/^A 2 - {E±i0+) 2 
and g%(E) = [g§{E) - gj(E)]f(E), where A is the su- 
perconducting order parameter (taken to be real). 

Now we have all the ingredients necessary to deter- 
mine the Green function in the dot and subsequently the 
spin and charge currents at both interfaces. The charge 
currents are given by 



h/r = 



16e 



dE Ti[T z i^ 



l/rl 



(6) 



Current conservation ensures that I = J; = — I r . 

The spin currents in the rotational frame are given by 



(7) 



In the laboratory frame, they decompose into a dc contri- 
bution along the precession axis, I a ,zi an d ac components 
in the perpendicular plane, I a ,x/y(t) — Ia,x/y cos ^ T 
I ay i x sin fit. Contrarily to the charge current, the spin 
currents do not need to be conserved. Eq. @ yields 



Ii + I r --^- 5 J dE Tr[(z x <r)T z g K ] 



= 0. 



(8) 



Thus, only the dc spin current along z is conserved. 

While our main interest are the FMR-induced cur- 
rents in the subgap regime of an F-dot-S junction, we 
first study the simpler case of an F-dot-N junction for 
comparison. For better readibility, in the following, we 
will normalize conductances by Gs = G; + G r and en- 
ergies by the Thouless energy E g = G^S/ (4ttGq). In 
particular, we introduce the dimensionless conductances 
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7a = Gq/Ge (a = l,r,m) as well as the dimensionless 
energies e = E/E g and uj = Q/(2E g ). 

A normal lead is described by setting A = in the 
above equations for g r . In the absence of superconduc- 
tivity, the retarded and advanced Green functions in the 
dot are trivial, g R ( A > = ±r z . The Keldysh component is 
obtained with the help of Eqs. ©, ©, and ©. There is 
a remarkable relation between the spin current at the left 
contact with the F lead and the charge current to lowest 
order in j m , 

I = m(t) ■ Ii(t) = [Ii tX sm6 + Ii tZ cos 8), 



ll 



ll 



(9) 

namely the charge current is proportional to the pro- 
jection of the spin current onto the instantaneous mag- 
netization axis of the barrier due to the spin-dependent 
conductance G m . That is, the charge current originates 
from two effects: (i) the rectification of the ac in-plane 
spin current pumped from the ferromagnet, and (ii) the 
conversion of the dc spin current along the z-axis into 
a charge current. It turns out that the two effects have 
opposite sign, and that the former dominates over the 
latter. Namely, we find 



GlE g Uj(j r + w 2 ) . 

L x — 7t- s — sin 8 cos 

2e 2 1 + lj 2 

J GlEg . 2 

Il,z = — ^2^ w 7rSin 8. 



(10a) 
(10b) 



Note that, in the limit 7 r < 71, the spin current along 
the z-axis is negligible. By contrast, in the limit 72 <C 
7 r , the two components are of comparable magnitude 
and almost complete cancellation between the competing 
effects takes place. 

The charge current reads 



I 



G m E g 7r7iW 3 sin 8 cos 8 
e 1+uj 2 



(11) 



At large precession frequency, lj ^> 1, the 
current scales linearly with frequency, / ~ 
(GiG r G m /2eG 2 ^)rism 2 8 cos8. In particular, in an 
open-circuit geometry, this would correspond to an 
FMR-induced dc voltage eV = (G m /2G^)Vt sin 2 8 cos 8 
in accordance with Refs. 0-11] ■ At w -C 1, spin- 
relaxation mechanisms induced by the tunnel coupling 
of the dot to the leads tend to suppress the effect. 

We now turn to the F-dot-S junction. In the subgap 
regime, the spin current at the interface with the super- 
conductor vanishes, I r =0. Thus, I z j = 0. However, an 
ac spin current is present at the interface with the fer- 
romagnet. Then, the Andreev charge current originates 
entirely from the rectification of this ac spin current. 

Restricting ourselves to energy scales much smaller 
than A, the Green function in the superconducting lead 



takes the simple form g. 



R(A) 



t x and g^ — 0. Taking 



7 m as a small parameter, we search for a perturbative 
solution of equation ([2]) in the form g = go + j m gi + • ■ • 
Due to the proximity effect, now the retarded and ad- 
vanced Green functions of the dot are modified as well. 
To zeroth order in j m , an explicit solution is given by 



*R(A) _ lrT x + [-i(e + cua z ) ± 7;]t z 
9o — 



^j 2 - (e + uja z ±i-fi) 2 



(12) 



Here, "f r is the effective minigap due to the coupling with 
the S lead [l6| , u an effective exchange field, and ji yields 
a broadening of the energy levels due to the coupling with 
the F lead. The Keldysh Green function can be cast in 
the form g$ = g R (p — (pg 1 



%A with 



<P = ./+ + /- cos* 



ll sin 8 

—5— oillVx-UVy) + COS 80 z 



(13) 



The function ip can be interpreted as a matrix- 
distribution function. Note that it does not depend on 
7 r as subgap electrons only thermalize with the F lead. 

To first order in j m , a solution which satisfies the nor- 
malization condition, g 2 = 1, is obtained in the form 
gi = c/qX — Xgo- For the advanced and retarded compo- 
nents, one finds X R ^ = ^f(sin8/2oj)[i^/ r /(e ± ili)r x + 
T z ](7 y . The Keldysh component can be decomposed as 
X K = X^t x + X*t Zi where Xf and Xf solve the cou- 
pled equations 



2>nX?-iu[a„X?] 



2sin#[cos# sin# /_ + f-\-(cr x ^yjh 

UJ 

2eX? -2iry r X? + u{*„X?} 
7i7 r sin " 



(14a) 



(14b) 



2i 



w( 7 2 + e 2 ) 



lif+(Jy - e/_ cos $(cos 8(7 x - smdcr z )]. 



Evaluating the current at the right interface, Eq. ([6|) 
yields I = — ry m G r /(16e) J dE Trf-r^gf]. Inserting 
the solution for iff and using the property g R (—e) = 
—&xT z go(e)a x T z , we obtain the current 



I = 



1 ifu 

2 7 2 +^ 2 



dc 



ef- 



il? + e 2 ){e + uj) 



(15) 



± 



-ji(e + uj) ± i(~ff - euj) 

Vlr ~ (e + W±i7;) 2 



where 1$ = (G m E g /e) sin 2 8 cos 8. The current as a func- 
tion of frequency for different values of 7; = 1 — j r is 
shown in Fig. [2] Simple analytic expressions can be found 
in different asymptotic regimes. 

In particular, at temperature T = and low frequency 
\cj\ <C 7, where 7 = (7 2 + 7 2 ) 1 / 2 , the FMR-induced cur- 
rent is given as 



10 r lhl,5 

3 7' 



(16) 
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n/(2E f ) n/(2E g ) 



FIG. 2: Andreev current induced by ferromagnetic resonance 
as a function of precession frequency for different values of 
7i = Gi/{G r + Gi). Here I = {G m E a /e) sin 2 <9cos#. 



FIG. 3: Induced current in the superconducting (dotted 
line) and normal state for 7; = 0.4. The thin line shows the 
contribution to the normal state current due to rectification 
only. 



The large power w 5 indicates the strong suppression of 
the effect. 

At large frequencies |w| 7, the current saturates. 
The frequency- independent value is given by 



-I sign(cj) x 



It 
1, 



Ir < Jl, 



(17) 



The saturation can be understood as Andreev processes 
become inefficient at energies larger than the minigap 

Depending whether the dot is more strongly coupled to 
the ferromagnet or to the superconductor, the crossover 
between these asymptotic regimes is different. If the dot 
is weakly coupled to the superconductor, 7 r <C 72, a 
smooth crossover happens at w ~ 1 > 7 r with a typi- 
cal current I/Iq ~ 7 2 . By contrast, if the dot is weakly 
coupled to the ferromagnet, 7; <C 7r , the crossover in the 
region u ~ 1 /2 is described by 



7 ;/(2v^fc), -l«^«- 7i , 
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where <5w = u — 1/2, with a typical current J/Jq ~ y/^fi 
at w = 1/2. 

While in the asymptotic regimes of cj very small or 
very large, the current is suppressed as compared to the 
normal state, there is in fact a wide intermediate regime 
where it may be strongly enhanced. Comparing Eqs. (fTTj) 
and JT7J), one notices that, if the dot is strongly coupled 
to the superconductor, in the regime 1/2 < 10 < r )^ 1 
the induced current in the superconducting state exceeds 
the induced current in the normal state. This effect may 
be understood due to the absence of a dc spin current 
along along the z-axis which leads to a strong suppres- 
sion of the effect in the normal state. Fig. [3] shows the 
current in the superconducting and normal state as well 
the contribution due to rectification only in the normal 
state. The ratio between the current in the supercon- 
ducting state and the latter contribution in the normal 



state reflects the ratio between Andreev and normal state 
conductances in an N-dot-S junction. 

So far we assumed that the magnetization in the ferro- 
magnet is uniform. However, boundary effects may lead 
to a suppression of the magnetization in the vicinity of 
the F/N interface. This would result in a different reso- 
nance frequency at the barrier than in the ferromagnetic 
reservoir and, consequently, in a tilt angle 9 b =/= 9 at 
FMR. The effect can be accounted for by replacing m 
with tub = (sm9s, 0, cos9b) in Eq. ([3]). In particular, 
at 9b = 0, the spin dependent conductance G m refers to 
the constant axis z. 

In the normal case, the relation between spin 
and charge currents, Eq. ([9]), now reads I = 
(2e7 TO 7 r /7; )rafl(i) • Ii(t). While at 9b = 9 the recti- 
fication of the in-plane ac spin currents always domi- 
nates over the conversion of the dc spin current along 
z into a charge current, this effect is completely sup- 
pressed at 9b = 0. As a consequence, at 9b = 0, the 
charge current, I = — (G m E g /e)ujj^ sin 2 9, has the op- 
posite sign compared to Eq. (fTTj) . In general, both ef- 
fects are important. The sign reversal occurs at tan#e = 
[1 - 7/w 2 /(7 r + us 2 )] tan 9. 

In the superconducting case, the dc spin current along 
z is always zero, and the charge current is due entirely 
to the rectification of the in-plane ac spin currents. As a 
consequence, we find that the charge current vanishes at 
9b = 0. The general result is obtained from Eq. (fT5|) by 
replacing I with Iq — (G m E g /e) sin#e sin 9 cos 9. 

In summary, we demonstrate that a subgap charge cur- 
rent in an F/S junction may be induced by ferromagnetic 
resonance. The effect is due to the rectification of ac 
spin currents generated by the precessing magnetization 
in the ferromagnet. In the normal case, a competing 
effect of conversion of a dc spin current into a charge 
current exists. This effect is absent in an F/S junction as 
the superconductor cannot carry a subgap spin current. 
As a consequence, the induced current in the supercon- 
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ducting state may be strongly enhanced as compared to 
the normal state. Interesting non-equilibrium phenom- 
ena should be expected in ferromagnetic Josephson junc- 
tions under ferromagnetic resonance conditions. 
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